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ABSTRACT 

Spitzer spectroscopy has revealed that ~ 80% of submm galaxies (SMGs) are starburst (SB) dom- 
inated in the mid-infrared. Here we focus on the remaining ~ 20% that show signs of harboring 
powerful active galactic nuclei (AGN). We have obtained Spitzer-IKS spectroscopy of a sample of 
eight SMGs which are candidates for harboring powerful AGN on the basis of IRAC color-selection 
(<S'8/jm/5'4.5pni > 2; 1.6. likely power-law mid-infrared SEDs). SMGs with an AGN dominating (> 50%) 
their mid-infrared emission could represent the 'missing link' sources in an evolutionary sequence in- 
volving a major merger. First of all, we detect PAH features in all of the SMGs, indicating redshifts 
from 2.5-3.4, demonstrating the power of the mid-infrared to determine redshifts for these optically 
faint dusty galaxies. Secondly, we see signs of both star- formation (from the PAH features) and AGN 
activity (from continuum emission) in our sample: 62% of the sample are AGN-dominated in the 
mid- infrared with a median AGN content of 56%, compared with < 30% on average for typical SMGs, 
revealing that our IRAC color selection has successfully singled out sources with proportionately more 
AGN emission than typical SB-dominated SMGs. However, we find that only about 10% of these AGN 
dominate the bolometric emission of the SMG when the results are extrapolated to longer infrared 
wavelengths, implying that AGN are not a significant power source to the SMG population overall, 
even when there is evidence in the mid-infrared for substantial AGN activity. When existing samples 
of mid-infrared AGN-dominated SMGs are considered, we find that S'8/im/'S'4.5^in > 1-65 works well 
at selecting mid-infrared energetically dominant AGN in SMGs, implying a duty cycle of ~ 15% if 
all SMGs go through a subsequent mid-infrared AGN-dominated phase in the proposed evolutionary 
sequence. 

Subject headings: galaxies: evolution — galaxies: high-redshift - galaxies: starburst — galaxies: active 
— infrared: galaxies — submillimeter 
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1. INTRODUCTION 

The era some 3Gyrs after the Big Bang {z ^ 2-2.5) 
coincides with a peak in activity in two important popu- 
lations: QSOs, which represent accretion onto supermas- 
sive black holes (SMBHs); and a population of extremely 
luminous, but highly obscured galaxies (e.g. Chapman et 
al. 2005; Wall, Pope & Scott 2008). The bulk of the lu- 
minosity of these obscured galaxies is emitted in the rest- 
frame far-infrared waveband. As a result of redshifting, 
they are most directly selected through their emission in 
the submillimeter (submm) or mm wavebands, typically 
in the atmospheric windows around 850 or 1100 /im, and 
so are termed submm galaxies (SMGs). The infrared 
luminosities (iin) inferred for SMGs from their submm 
emission are highly uncertain, but assuming, as appears 
to be the case, that they follow the far-infrared-radio cor- 
relation for local starburst (SB) galaxies, then their typ- 
ical luminosities will be Lir ~ 10^^-10^'^ L0 (e.g. Kovacs 
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et al. 2006; Murphy et al. 2009). Thus this population 
may contain some of the most luminous galaxies in the 
Universe comparable in luminosity to QSOs. 

The increasing availability of precise redshifts for sam- 
ples of SMGs (e.g. Chapman et al. 2005; Eales et al. 

2009) has allowed their properties to be studied in de- 
tail. SMGs are strongly clustered (Blain ct al. 2004), 
massive (M* > 10" Mq, Borys et al. 2005; Hainline 
et al. 2009), and gas rich (/gas ~ 0.3, Prayer et al. 
1998; Greve et al. 2005; Tacconi et al. 2006, 2008) 
systems which are known to harbor (apparently) low- 
luminosity Compton-thin AGN (i.e. A^h < 10^^ cm^^, 
Takata ct al. 2006; Alexander ct al. 2005a,b, 2008a) and 
(apparently) strong star-formation (SF) activity (SFR ~ 
lOOOMgyr-i; Swinbank ct al. 2004; Chapman et al. 
2005). Many of these properties, and the similarity 
between the redshift distributions of QSOs and SMGs 
(Chapman et al. 2005), support a link between SMGs, 
QSOs and the formation phase of massive elliptical galax- 
ies (e.g. Lilly et al. 1999; Archibald et al. 2002; Stevens et 
al. 2005; Coppin et al. 2008b). In the high-redshift inter- 
pretation of the evolutionary sequence first presented by 
Sanders et al. (1988), SMGs would trace an infrared ul- 
traluminous phase followed by a short 'transition phase' 
where the galaxy would display a mix of SF and obscured 
AGN activity before evolving into an optically luminous 
QSO. Studying the relative SF and AGN activity in these 
'transition' or 'missing link' sources and performing a 
comparison to typical SF SMGs can thus provide impor- 
tant insight on the validity of this evolutionary sequence. 
Here we focus on the comparison of the mid-infrared 
spectral properties between existing samples of typical 
SMGs and a new sample of candidate 'transition phase' 
SMGs. The benefit of undertaking this energy audit in 
the mid-infrared is that the spectra can provide measure- 
ments of both SF and AGN activity simultaneously (see 
Pope et al. 2008a and Menendez-Delmestre ct al. 2009; 
hereafter P08 and MD09), providing a complementary 
approach to radio (e.g. Ibar et al. 2010) and X-ray stud- 
ies of SMGs (e.g. Alexander et al. 2005b; Laird et al. 

2010) . 

The Spitzer InfraRed Spectrograph (IRS; Houck ct al. 
2004) era enabled the study of the energetics of « 45 
SMGs with 24 jim flux densities as low as ~ 0.1 mJy and 
with redshifts as high as z = 2.6, demonstrating that 
accurate redshifts for dust-enshrouded (and sometimes 
optically invisible) galaxies can be obtained (e.g. Lutz et 
al. 2005; Valiante et al. 2007; Menendcz-Delmcstre ct al. 
2007; MD09; P08). The majority of the IRS work has 
confirmed that SMGs are primarily SB-dominated sys- 
tems, with hot dust continuum from an AGN contribut- 
ing at most 30% of the mid-infrared luminosity. Only 
about 15% of blank-field SMGs appear to be continuum- 
dominated in the mid-infrared (i.e. > 50% contribu- 
tion; P08), with InfraRed Array Camera (IRAC; Fazio 
et al. 2004) 8.0 /xm to 4.5 /im color ratios of Ss/S4,5 > 
2. Mid-infrared continuum-dominated SMGs are poten- 
tially an important sub-population of SMGs representing 
the 'missing link' sources in the proposed evolutionary se- 
quence of Sanders et al. (1988), but being the minority of 
this luminous population, have not yet been studied in a 
systematic or statistically robust way. Thus, an impor- 
tant question to address is: Are the energetics of these 
composite SB/ AGN objects with submm emission dom- 



inated by SBs or by AGN, even when the IRAG colors 
indicate that an AGN is likely present? 

Other relevant IRS samples include near-infrared- 
selected SBs, X-ray-selected AGN (Weedman et al. 
2006), and Spitzer-selected 0.3 < 2 < 3 ULIRGs (Yan et 
al. 2007; Sajina et al. 2007; Farrah et al. 2008; Dasyra et 
al. 2009; Hernan-Caballero ct al. 2009). AGN-dominated 
sources in these samples tend to reside in a distinct pa- 
rameter space in Spitzer color-color diagrams, which is 
consistent with color-rcdshift evolution tracks of well- 
known local AGN (sec e.g. Fig. 1 ; Ivison et al. 2004; 
Ashby et al. 2006; Hainhnc ct al. 2009). The reason why 
Ss/S4,5 > 2 should locate mid-infrared AGN-dominated 
sources is simple: seeing an enhanced 8 /im flux density 
compared with 4.5 /im is expected if a source has sig- 
nificant thermal power-law emission from an AGN ac- 
cretion disk, which can dilute both the PAHs and thc^ 
H-opacity minimum (1.6 /um stellar bump, which these 
channels trace). Although some contamination is ex- 
pected: SF-dominated sources can also show enhanced 
'S'8/'5'4.5 at z > 4 (when these channels begin to sample 
over the peak of stellar photospheric emission) aswcll as 
at 2; ^ 1 (since these channels are not yet climbing up 
the restframe 1.6 /um stellar bump). The color cut is thus 
appropriate for separating mid-infrared SB- and AGN- 
dominated SMGs from z ~ 1-4, which is well matched 
to the known redshift distribution of the SMG popula- 
tion. 

We have thus selected eight SMGs from the Submil- 
limeter Common-User Bolometer Array (SCUBA; Hol- 
land et al. 1999) HAlf Degree Extragalactic Survey 
(SHADES; Morticr et al. 2005; Coppin et al. 2006; 
Austermann et al. 2010) within this relatively imexplored 
Ss/S4,5 > 2 parameter space that are likely harboring 
AGN. These were targeted with the IRS in order to deter- 
mine the relative contribution of power-law/ AGN emis- 
sion versus PAH/SF emission to their power output, en- 
abling a comparison to similar IRS samples of more typ- 
ical SF-dominated SMGs. We use these data to obtain 
independent redshift estimates as well as to test if the 
IRAC-color criterion is a secure means of pre-selecting 
SMG counterparts with an enhanced AGN component 
compared to typical SMGs - the 'missing link' sources 
in the proposed evolutionary sequence we wish to inves- 
tigate. 

This paper is organized as follows. The sample selec- 
tion, Spitzer-IRS observations, data reduction, and anal- 
ysis approach arc described in § 2. In § 3 wc present the 
main results of the IRS spectroscopy, inchiding redshifts, 
spectral decomposition and AGN classification, and full 
SED fits to determine their total infrared luminosities. 
We discuss the implications that our results and those 
of other IRS SMG studies have on the role that SMGs 
play in galaxy evolution in the framework of the proposed 
evolutionary sequence in § 4. Finally, our conclusions are 
given in § 5. 

Here, we discuss the observed 'mid-infrared' spectral 
properties of our SMGs probed by the IRS which, at 
our source redshifts of z > 2.5, traces ^ 4 10 /im in the 
rest-frame. All magnitudes in this paper are on the AB 
system, unless otherwise stated. We adopt cosmological 
parameters from the WMAP fits (Spergel et al. 2003): 
= 0.73, = 0.27, and Hq = 71 kms"! Mpc"!. 
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Fig. 1. — Spitzer color-color diagram as an AGN diagnostic (see e.g. Ivison et al. 2004; P08). The dotted and dot-dashed curves show 
the positions of Mrk231 (an AGN) and M82 (a SB), respectively, as a function of redshift (redshift is indicated by the numbers along the 
tracks). The smallest filled and open circles represent robust (P < 0.05) and tentative (P > 0.05) 24 ^m SMG counterparts, respectively, 
for our parent sample of blank- field SHADES SMGs (Ivison et al. 2007). For comparison, we have also plotted the GOODS-N SMG sample 
(blue squares; Pope et al. 2006), and also the Chapman et al. (2005) radio-selected spectroscopically-confirmed sample of SMGs (green 
triangles). SMGs from these parent samples that have been observed with the IRS (P08; MD09; this work) are indicated by large versions 
of the corresponding symbols. Of these IRS-observed SMGs, those with > 50% contribution from continuum (AGN) emission to the 
mid-infrared luminosity are circumscribed and tend to lie far from the M82 SB sequence. Although our sample was initially selected on the 
basis of 5g/54.5 > 2 (thin dashed vertical line), when all the results from the literature are taken into account a selection of SS/S4.5 > 1.65 
thus seems to pick out AGN fairly efficiently, with some small amount of scatter over this boundary (thick dashed vertical line; see § 4). 
Extrapolating from this result implies that 15% of blank-field SMGs will likely have an enhanced dominant AGN component in the 
mid-infrared. 



2. SAMPLE SELECTION, OBSERVATIONS, DATA 
REDUCTION, AND ANALYSIS 

2.1. Sample Selection and Multi-wavelength Properties 

SHADES provides a suitably large parent sam- 
ple of SMGs from which we can select a repre- 
sentative and large enough sample of SMG AGN- 
candidates. SCUBA surveyed ~ 0.25 dcg^ at 850 /im 
to an rms of ~ 2mJy, uncovering 120 SMGs > 3.5 cr 
in the Lockman Hole East (LH) and the equatorial 
Subaru- XMM Deep Field (SXDF), centred at J2000 
RA=2i^l7'"57!5, Dec=-5°00'18'.'5 and RA=10^52™26!7, 
Dec=57°24'12'.'6, respectively (Coppin et al. 2006). The 
SHADES fields were subsequently mapped to their full 
~ 0.5 deg^ coverage by the Astronomical Thermal Emis- 
sion Camera (AzTEC; Wilson et al. 2008) at 1100 /im 
to an rms depth of ~ 1-1.5 mJy, yielding 114 SMGs 
above 3.6 ct (Austermann et al. 2010), a population akin 
to the well-known submillimeter or SCUBA galaxies 
(e.g. Chapin et al. 2009). A detailed comparison of the 
overlap in the SCUBA and AzTEC SHADES catalogs 
will be presented in Negrello et al. (in preparation). 

Ivison et al. (2007), Ibar et al. (2009, 2010) and Aru- 
mugam et al. (in preparation) have identified 1.4 GHz 
Very Large Array (VLA), 610 MHz Giant Metre- wave 
Radio Telescope (GMRT), and/or Mid-Infrared Pho- 
tometer for Spitzer (MIPS; Rieke et al. 2004) 24 ^m 
counterparts for 65% of the SHADES SMGs. Here we se- 
lect SMGs with statistically robust 24 /im counterparts 
with > 0.2 mJy (ensuring that we achieve ade- 

quate signal-to-noise, S/N, on the continuum in the IRS 
spectra) and observed IRAC colors of Ss/S4,5 > 2 (Dye 
et al. 2008; Clements et al. 2008; see Fig. 1). Our final 



sample is comprised of eight sources: four SMGs from 
each of the SHADES SCUBA and AzTEC surveys. 

Our selection criteria apply to approximately 15-25%^ 
of the SHADES SMGs (or equivalently to 7-13% of 
SMGs with 24 /im counterparts). Note that including 
the 'non-robust' 24 /im IDs in these estimates gives sim- 
ilar fractions. In addition, deeper radio and 24 /im data 
would likely yield identifications for all of the SMGs; and 
there is no evidence in support of SMGs fainter than 
200 /i Jy at 24 /im being fundamentally different than per- 
haps lying at slightly higher redshifts (e.g. Pope et al. 
2006). Thus, the IRAC color-selected subset of SMGs 
studied here should be representative of ^ 15% of the 
blank-field SMG population. We discuss the possibility 
of more subtle selection effects in § 4.1. 

We also observed a 24 /im counterpart (LOCK850.41- 
2) of an interesting SMG with multiple robust radio and 
MIPS counterparts, where the other MIPS counterpart 
(LOCK850.41-1) was previously observed by MD09. We 
present these data separately in Appendix A, since this 
source does not make the original color-cut criterion for 
our sample. 

We make use of additional multi-wavelength data-sets 
available for the 24 /im SHADES SMG counterparts in 
order to further interpret our results in the wider context 
of the entire SMG population and other high-redshift 
galaxy populations. Available optical-to-radio photom- 
etry for each source is given in Table 2 and was retrieved 
from Dye et al. (2008), Furusawa et al. (2008), and 

^ The range in this fraction merely represents the fact that sev- 
eral of the SMGs are undetected 8 ^m and so they could lie on 
either side of the IRAC color-cut boundary in Fig. 1. 
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Morticr ct al. (in preparation; Subaru BRi' z), Lawrence 
et al. (2007) and Warren et al. (in preparation; J and K 
from the DR3 UKIRT Infrared Deep Sky Survey release 
- UKIDSS, Egami et al. (priv. comm.), and Dunlop et 
al. (in preparation; 3.6-8 /xm Spitzer-IRAC) . 

The SHADES fields also possess moderately sensitive 
XMM-Newton imaging with exposures of 673 ks in the 
LH (Brunner et al. 2008), and 18-83 ks in the SXDF 
(Ueda et al. 2008). In the deeper X-ray data, even X-ray 
non-detections, when combined with IRS data, can be 
sufficient to imply 1-2 orders of magnitude of obscuration 
at rest-frame 2- lOkeV and can identify Compton-thick 
AGN. Given the combined positional uncertainties of the 
24 /xm and X-ray sources (yielding a combined positional 
uncertainty of 3 4"), we search for X-ray counterparts 
to the SMGs using a 4" search radius. For any X-ray 
detections, wc extract a 0.5-2 keV flux and calculate the 
rest-frame 2-10 keV X-ray luminosity, assuming a spec- 
tral slope of r = 1.4 (the average X-ray spectral slope 
for sources of this approximate X-ray fiux; see Fig. 8 
of Alexander et al. 2001) to make a small K-correction 
(;C 30%). For the X-ray non-detections, we have derived 
3 a upper limits. The X-ray fluxes and flux upper limits 
are given in Table 3. 

It happens that spectroscopic redshifts have not been 
obtained for the majority of the SMG sample, undoubt- 
edly due to the faintness of the corresponding opti- 
cal counterparts (sec Table 2). All of our SMGs are 
i?AB > 24 (except for AzLOCK.lO which is i?AB = 22.65, 
and which has not yet been attempted in a spectroscopic 
follow-up program, to our knowledge), for which Chap- 
man et al. (2005) flnd an increased spectroscopic rcdshift 
failure rate. Therefore it is not a surprise that only one 
of our targets has a published tentative optical spectro- 
scopic redshift available for comparison with the IRS- 
derived redshift (see Appendix A). 

2.2. Spitzer-IRS Observations 

The Spitzer-IRS observations (PID 50183) were taken 
in spectral mapping mode, with the target being placed 
at six positions (separated by 20") along the slit using 
multiple cycles of the longest ramp setting of 120 s to 
maximize S/N in unit time, while ensuring internal ro- 
bustness against cosmic rays and the identification of 
rogue pixels (see Teplitz et al. 2007). We observed the 
SMGs using Long-Low 1 (LLl; 19.5-38.0 /xm) and Long- 
Low 2 (LL2; 14.0-21.3/xm) to ensure that the spectra 
cover a significant fraction of the full range of PAHs (from 
6-17 /xm in the rest-frame). Integration times were cal- 
culated on a source-by-source basis using the observed 
24 /xm fluxes to yield a S/N of > 3-4 in the contin- 
uum near the band centres (i.e. similar-quality spectra 
to those in P08 and Mencndez-Delmestre et al. 2007, 
2009). We used high-accuracy blue IRS peak- up acqui- 
sition on nearby isolated moderately bright Two Micron 
All Sky Survey stars (2MASS; Skrutskie et al. 2006). To 
verify the calibration in the LL2 spectra we also obtained 
16 fim IRS peak-up imaging with 2 x 30 s cycles in a 5- 
point dither for each target for an additional 1.7 hrs in 
total. Table 1 gives our target list and integration times 
for the spectroscopy. Our observations were obtained 
in 2008 May, June (LH targets), and September (SXDF 
targets) in a total of 58.4 hrs. 



2.3. Data Reduction 

We begin with the Spitzer pipeline (Version S18) ba- 
sic calibrated data. In reducing the IRS data, we follow 
the same approach used to reduce IRS spectra of faint 
sources outlined in P08. To summarize, we identify and 
clean rogue pixels in the 2D files, fit and subtract latent 
build-up on the arrays, perform sky subtraction by cre- 
ating a normalized 'supersky' for each AOR in which all 
other sources have been masked, and co-add all the 2D 
sky-subtracted data files for each nod position. We ex- 
tract the ID spectrum from the 2D co-added files for each 
nod using the the optimal-extraction mode in SPICE^. 
Along with the target spectrum we also extract a residual 
sky spectrum from each 2D co-add file, which we use to 
determine uncertainties on our final ID spectra. We trim 
all the spectra to < 35.0 /zm in order to restrict the anal- 
ysis to the reliable and least noisy portion of the data. 
The observed spectra are shown in Fig. 4. 

We process the IRS 16 /im peak-up images (PUIs) with 
the SSC pipeline S18. We then correct for latent charge 
accumulation in the PUIs by subtracting the mode flux 
of the central 30x45 pixels from each individual expo- 
sure. We eliminate residual sky background by sub- 
tracting from each science exposure a median sky, cre- 
ated from median-collapsing all charge-corrected expo- 
sures and scaled to match off-source background level 
for each exposure. All sky-subtracted, charge-corrected 
science exposures are mosaicked into a single image us- 
ing the standard MOPEX mosaicking pipeline, keeping 
the native PUI resolution scale of 1.8"/pix. We derive 
16 /im photometry by performing point-response func- 
tion (PRF) fitting analysis with the MOPEX/ APEX 
single- frame pipeline, relying on the blue PUI PRF^. We 
use 'peak' thresholding as part of the image segmenta- 
tion algorithm, setting the detection threshold to a 2a 
level within a 5x5 pixel fitting area to ensure adequate 
detection of faint sources. To characterize the uncertain- 
ties in the derived integrated fluxes we rely on the inte- 
grated noise within an on-source 5x5 pixel aperture in 
the uncertainty mosaic. The 16 /im imaging photometry 
is reported in Table 2 and shown in Fig. 2. 

We verify the calibration of our final ID spectra by 
comparing to the MIPS 24 /xm and IRS/PUI 16 /xm 
fluxes. In all cases the spectra are consistent with the 
imaging photometry within the uncertainties. 

2.4. Analysis 

Our aim is to determine if our mid-infrared spectra 
show signs of AGN-dominance, and if so, to then compare 

these spectral properties with the overall SMG popula- 
tion. We have thus followed a similar analysis approach 
to P08 in order to facilitate a direct unbiased comparison 
with existing Spitzer-IRS samples of SMGs. We have or- 
ganised the paper so that the following subsections con- 
tain all the salient details of our analysis approach, with 
the results presented separately in § 3. In summary, in 
§ 2.4.1, we identify the PAH features present in our spec- 
tra and determine redshifts. We then decompose each 
individual spectrum into SF and AGN components in 
§ 2.4.2 in order to classify the fractional contribution 
of AGN to the mid-infrared luminosity. In § 2.4.3 we 

^ Available at http;//ssc. spitzer. caltech.edu/postbcd/spice.html 
^ Available at http://ssc.spitzer.caltech.edu/irs/puipsf 
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Fig. 2. — New 16 /jm 10" x 10" imaging (roughly the size of 
the aperture of the IRS observations) centered on the 24 /im po- 
sitions of the nine SMG targets (in order of increasing redshift; 
see Table 3; left-right/top-bottom), with VLA 2(t,3(t,... contours 
overlaid. Source names are given in each panel (see Table 1). The 
sample has been selected on the basis of being 24 fira detected and 
IRAC color-tuned in order to pick out potential AGN-dominated 
SMGs, although consequently the majority of the targets also ap- 
pear to be radio-bright, which is consistent with an AGN-biased 
sample. XMM-Newton detections are shown by the purple 'X' 
symbols, and only exist for the minority of the sample - deeper 
X-ray data would likely yield a higher detection fraction, unless 
the sources are Compton-thick (see text). 
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Fig. 3. — An example of the spectral line decomposition of an 
SMG IRS spectrum, used to compute the redshift for the source 
and to measure the line profile features (see § 2.4.1). The solid black 
curve is the raw (unsmoothed) IRS spectrum, while the shaded 
region represents the associated 1 a noise from the sky background. 
The red dashed curve shows the best-fit SED to the mid-infrared 
spectrum, comprising a strong continuum component (blue dot- 
dashed power-law) and PAH features (blue dot-dashed Lorentzian 
profiles) at 6.62, 7.71, and 8.61 fj,m (indicated by vertical lines). 

extrapolate our results to determine the full mid-to-far- 
infrared SEDs in order to estimate the AGN contribution 
to the bolomctric luminosity of these systems. 

2.4.1. Redshift Determination and Spectral Line 
Measurements 

In order to derive the source redshifts, we simul- 
taneously fit a power law continuum component and 
Lorentzian profiles to the spectra using a Levenberg- 
Marquardt least-squares method, including the residual 
sky errors in the fit determination (e.g. Sajina et al. 
2007), assuming the following rest- frame wavelengths for 
the main PAH features that are visible: 6.22, 7.71, and 
8.61 /im (Draine & Li 2007). We do not impose any joint 
constraints on the line ratios or centres (i.e. to all lie at 



the same redshift), as we aim to obtain individual red- 
shift estimates from each PAH feature. See Fig. 3 for an 
example of these fitting results. Our fitting routine re- 
turns the formal 1 a statistical error estimates for each set 
of best-fitting parameters derived from the square root 
of the diagonal elements of the covariance matrix, and 
we adopt these errors in the ensuing analysis. We verify 
that these formal errors are reasonable through Monte 
Carlo simulations similar to those described in P08. 

The final redshift for each SMG is calculated by tak- 
ing a weighted average of the individual redshift esti- 
mates from the visible PAH lines (although we neglect 
the 8.61 /im PAH feature in the redshift estimate be- 
cause it appears relatively weak and significantly blended 
with the much stronger neighboring 7.71 /im feature). 
The total error on the redshift is calculated from la = 
{J27=i ^pah )^^^'^ ^ith an additional 'centroiding error' 
(calculated as the maximum spread in the line centering 
from the Monte Carlo simulations) added in quadrature. 
For sources with only one PAH line, the 1 a redshift er- 
ror is taken to be the centroiding error derived from the 
simulations (the more conservative choice). Note that 
the redshift errors range from 0.01-0.05, consistent with 
the similar-quality spectra in P08 and MD09. For com- 
parison, we use the Brandl et al. (2006) SB template 
with fixed PAH linewidths and positions and find that 
this gives redshift estimates consistent with our method 
(which is essentially equivalent to performing a combined 
fit of all the lines simultaneously, but with a more accu- 
rate measure of the uncertainty of the redshift). 

Once redshifts have been obtained (see Table 3), we 
calculate the rest-frame PAH line equivalent widths 
(EWs) and luminosities from our Lorentzian model fits, 
propagating through the redshift and parameter errors 
from the fits (see Table 4). 

2.4.2. Determining the Relative AGN and SF 
Contributions in the Mid-Infrared 

We determine the relative AGN and SF contributions 
by fitting each spectrum with a model comprised of three 
main components (extinction, PAHs and a continuum 
power-law) of the form F^, = cq v~^^ e"^^"^" -|- C2 /j/,m82 
using a minimization approach (see P08)''. The AGN 
emission is characterized by a power-law with both the 
normalization and slope as free parameters, and T9.7 is 
obtained from the Draine (2003) extinction curves. We 
assume that all of the PAH line emission is powered by 
a SB (e.g. Rigopoulou et al. 1999; Laurent et al. 2000), 
and the SB /PAH emission is fitted using the mid-infrared 
spectrum of M82 from Forster Schreiber et al. (2003), 
which is known to fit and describe the mid-infrared spec- 
tra of the SB-dominated SMG population well (e.g. P08). 
We derive the fraction of the mid-infrared luminosity 
that comes from the continuum component by integrat- 
ing the continuum portion of the best-fitting model and 
dividing by the area under the total fit. The use of differ- 
ent templates was explored by Sajina et al. (2007) and 
P08 (e.g. using NGC7714 in place of M82), and these 
authors find that it makes little difference on the final 

In contrast to the more typical SMGs from P08, we found that 
the fitting never required additional extinction of the PAH compo- 
nent, so this term has been omitted from our model (cf. eqn. 1 in 
PCS). 
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Fig. 4. — IRS spectra of our AGN-candidate sample of SMGs, listed in order of increasing redsliift and shifted to rest-frame wavelengths 
in order to facilitate an easy inter-comparison of the spectral features. Source names are given in each panel (see Table 1). The solid 
black curves are the raw (unsmoothed) IRS spectra, v^hile the shaded region represents the associated la noise from the sky background 
(see § 2.3). The red dashed line is the best-fitting model (see § 2.4.2), which is composed of an M82 PAH template (green dot-dot-dot- 
dashed line) and a power-law component (blue dot-dashed line) with extinction applied. The continuum (AGN) fraction for each source 
is calculated from the spectral decomposition (see Table 3). The spectra and fitting reveal that our sample of candidate AGN-dominated 
SMGs all show PAH features atop a range of mid-infrared continuum slopes, indicating a mix of SB and AGN properties in the sample. 



outcome of the decomposition, especially given the small 
wavelength range covered by the spectra. In addition, 
we perform the decomposition using the Brandl et al. 
(2006) SB template in place of M82 and find that they 
yield comparable fractions. 

2.4.3. Full SED Fitting and Determination of Lir 

One of the main goals is to determine if our sample of 
SMG mid-to-far-infrared SEDs resemble that of typical 
SMGs or if an additional AGN component is required 
to fit the data. With this goal in mind, using the same 
simple templates and approach as P08 will enable a di- 
rect and fair comparison between the two samples. We 
have thus used the full mid-infrared-to-radio composite 
SED from P08, comprised of their mid- infrared SMG SB 
composite spectrum spliced with the best-fitting modi- 
fied Chary & Elbaz (2001) template (CEOl hereafter), to 
represent the SF component, and we have used the SED 
of Mrk231, consisting of its IRS spectrum (Armus et al. 
2007) spliced together with the CEOl template that best 
fits the far-infrared/submm photometry of Mrk231, to 
represent the AGN component. Although the resulting 
decomposition will thus depend heavily on the templates 
we have assumed, these choices will allow us to determine 
the relative differences between the spectral decomposi- 
tion of our sample of AGN-dominated SMGs and of more 
typical SB-dominated SMGs from P08. Nevertheless, we 
discuss and explore other choices for the SB and AGN 
templates § 3.4 to verify the dependence of our results 
on the assumed templates. 

We fit linear combinations of these templates to the 
data, including the mid-infrared IRS spectra, the 850 or 
1100 /xm photometry, and, where available, 350 /im pho- 



tometry from Coppin et al. (2008a) using a Levenberg- 
Marquardt least-squares approach. We also verify that 
available 70 and 450 /im 3 a upper limits from Hainline 
et al. (2009) and Coppin et al. (2006), respectively, are 
not violated in the final fits. We then integrate the total 
best-fitting linear combination of the SB and AGN SED 
components in the rest-frame from 8-1000 /im to obtain 
the total infrared luminosities, Xir, and calculate the 
relative fractions and luminosities of the SB (Lge) and 
AGN (Lagn) components to the total infrared luminos- 
ity (see Table 4). Note that the full SF-SMG composite 
SED from P08 already contains a 30% AGN contribution 
in the mid-infrared (from 5-11.5/im in the rest-frame), 
which is taken into account in our calculations of Xagn 
and LsB- 

After the level of AGN contamination in each SMG has 
been determined in this way, we calculate the SFRs for 
our objects from Lsb following Kennicutt (1998), assum- 
ing a starburst less than 100 Myr in age with a Salpeter 
(1955) initial mass function. 



3. RESULTS 

We now present the results of our experiment, includ- 
ing the source spectra and redshifts (§ 3.1), the spectral 
decomposition and a comparison with the X-ray view in 
order to gather a complete census of the energetics of the 
sample (§ 3.2), the construction of a composite AGN- 
dominated SMG (§ 3.3), and finally the full SED fits 
(§ 3.4). We stress that we have followed similar recipes 
to those of P08 in order to facilitate a direct comparison 
with samples of more typical SF-dominated SMGs (see 
§ 2.4). 
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3.1. Source Spectra and Redshifts 

We present the final reduced spectra and fits in Fig. 4. 

At a glance, the spectra reveal that our sample of SMGs 
possess PAH features atop a range of mid-infrared con- 
tinuum slopes. Despite the strong continuum in many of 
the spectra, measurable PAH features are apparent in all 
of our spectra, and we use these to extract redshifts (see 
§ 2.4.1 and Table 3). Near-infrared spectroscopic red- 
shifts have recently been obtained for SXDF850.02 and 
SXDF850.37 (Alaghband-Zadch et al. in prep.) of z = 
2.518±0.001 (showing narrow lines) and z = 2.769±0.001 
(showing broad Ha and Oiii5007), respectively, in agree- 
ment with the IRS-derived redshifts. Finally, we note 
the suggestion by Ivison et al. (2005) that a galaxy at 
^^spoc — 2.15, 3" west of LOCK850.1 may be associated 
with the SMG. Our Spitzer-IRS result suggests other- 
wise, and we discuss this source further in Appendix A. 

The redshift distribution of our color-selected sample 
of eight objects ranges from z = 2.5-3.4, with a median 
of 2.76. This is noticeably shifted to higher redshifts 
compared to the P08 and MD09 samples, which have 
a median redshift of 2.0, which is more consistent with 
the known redshift distribution for radio-identified SMGs 
with 5850 > 5mJy and z ~ 2.2 (Chapman et al. 2005). 
In § 4.1, we suggest that the cause for this overall red- 
shift difference between the sub-samples is likely due to 
selection effects. 

3.2. Spectral Decomposition of Individual Galaxies 

3.2.1. Relative AGN/SF Contribution in the 
Mid-Infrared 

Using the fitting procedure described in § 2.4.2 to de- 
termine the relative AGN/SF contribution in the mid- 
infrared, we have tabulated the continuum (AGN) frac- 
tion for each SMG in Table 3. We find a wide range of 
continuum fractions in our sample - from sub-dominant 
to significant - with a median continuum contribu- 
tion to the mid-infrared of « 56%. SXDF850.02 and 
SXDF850.37 have recently been spectroscopically con- 
firmed at restframe optical wavelengths by Alaghband- 
Zadeh et al. (in prep.), and the restframe optical line 
widths are in broad agreement with our mid-infrared 
AGN fractions of 10% (where narrow lines suggest neg- 
ligible AGN contribution) and 62% (where broad Ha 
and Oiii5007 lines suggest the presence of an AGN), re- 
spectively. 5/8 (62%) SA4Gs in our sample arc clearly 
continuum-dominated sources in the mid-infrared (> 
50% of the mid-infrared luminosity). For comparison, 
only 2/13 (15%) SMGs from P08 have > 50% AGN con- 
tribution in the mid-infrared, and the median AGN con- 
tribution of their sample (including those with upper lim- 
its) is « 34%. 

Although deeply-embedded Hii regions could poten- 
tially be contributing to the hot dust continuum that 
we here associate with an AGN contribution (making 
the fractions quoted above merely upper limits to the 
AGN contribiition) , we argiie below that the dust heated 
by stars is likely a minor contribution. In Fig. 5 we 
have plotted our 3.6-24 /xm photometry in the rest- frame 
to search for any evidence of the H-opacity minimum 
(1.6/xm stellar bump), as traced by the 4.5 and 5.8 /im 
channels given our source redshifts, which is normally 
visible in typical SF SMGs (e.g. Hainline et al. 2009; 




1 10 
Rest wavelength [\iva\ 



Fig. 5. — The rest-frame near-infrared SED of our sample of 
Ss/S4„5 > 2 SMGs (open cicles) compared with a. fourth-order 
polynomial fitted to 58/54.5 < 2 (SB-dominated) SMGs from Mes- 
senger et al. (in preparation). This simple comparison shows that 
on average our sample of S8/S4.5 > 2 SMGs are dominated by 
power-law continuum emission which swamps the stellar emission 
bump normally seen in typical SB-dominated SMGs. 

Messenger et al. in preparation). The rest- frame broad- 
band photometry of our sample of SMGs reveals a strong 
power-law shape from 1-3 /xm. This exercise reveals that 
on average our sample of Ss/S^.s > 2 SMGs shows no 
1.6 /zm stellar bump, whereas more typical SB-dominated 
SMGs (58/S'4.5 < 2) do (sec also Yun et al. 2008). We 
thus interpret the Sa/Si,^ > 2 observed color excess as 
due to hot dust emission (500-1000 K) in the immediate 
vicinity of an active nucleus (Laurent et al. 2000). By ex- 
tension, this AGN emission would also dominate at the 
longer wavelengths probed by the IRS spectra (< 10 /im), 
causing the power-law shape that we observe. Tentative 
evidence for AGN emission in these SMGs is apparent at 
radio wavelengths in the form of enhanced 1.4 GHz emis- 
sion compared to 610 MHz data (see § 3.4.1). We can 
directly test how likely it is that the AGN is producing 
the steep continuum emission in the mid-infrared spec- 
tra by examining the combination of the mid-infrared 
and X-ray data, as we now describe. 

Using the mid-infrared continuum fits from § 2.4.2, we 
calculate rest-frame 6 /im luminosities due to the AGN 
component and compare these with the X-ray luminosi- 
ties (Table 3). Comparing the X-ray and mid-infrared 
data can potentially tell us about the intrinsic luminos- 
ity of the AGN and can reveal whether any of the sources 
are Compton-thick by comparing the estimated column 
densities to other samples of SMGs and local AGN. Fol- 
lowing Alexander et al. (2008b), in Fig. 6 we have plotted 
the rest-frame 2-10 keV versus 6/im luminosities for our 
sample compared to AGN in the local Universe and to 
other SMGs. Assuming the X-ray -to-mid- infrared lu- 
minosity relationship for local AGN (Lutz et al. 2004) 
holds at high-redshift, the combination of the intrinsic 
X-ray luminosities and 6 /xm luminosities of some of our 
objects suggests that they may be Compton-thick (i.e. 4 
have upper limits close to the lower dotted line in Fig. 6). 
While the X-ray data are not sufficiently deep (cf. 2 Ms 
Chandra data in GOODS-N, e.g. Alexander et al. 2008b) 
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Fig. 6. — Rest-frame 2-10 keV luminosity derived from the X-ray 
versus 6 luminosity derived from the mid-infrared power-law 
component fit (filled circles, where arrows indicate X-ray 3 a upper 
limits). For comparison, we overplot Spitzer-IRS z ~ 2 AGN with 
optical spectroscopy from Alexander et al. (2008b) (open squares) 
and the SMG CI from P08 (open circle) which have properties 
suggesting they are Compton-thick . We have also overplotted the 
intrinsic (absorption corrected) X-ray-to-mid-infrared luminosity 
relationship for local AGN (Lutz et al. 2004; the shaded region) 
and, assuming that it holds at high-redshift, the combination of 
the observed (not corrected for absorption) X-ray luminosities and 
6 ^m luminosities of some of our objects suggests that they could 
be Compton-thick (see Fig. 4 of Alexander et al. 2008b; the lower 
dotted line refers to the observed luminosity ratio expected for a 
typical AGN that is absorbed by Afjj — 10^* cm~^ using the model 
presented in Alexander et al. 2005b). Note that CI from P08 is 
also continuum-dominated in the mid-infrared, but it possesses a 
much lower Lx upper limit since it lies in a much deeper X-ray 
field. 

to conclusively determine whether or not the AGN are 
Compton-thick, the mid-infrared spectra show that AGN 
are dominating in these sources (in the form of strong 
continuum emission), and the non-detections are consis- 
tent with sources with being absorbed at X-ray energies 
(i.e. TVh ^ 10^2 cm-2). Indeed since 7/8 SMGs are also 
optically faint with R > 23, the majority of the sam- 
ple appear to be obscured to some degree. Compared to 
other SMGs, most of our objects are consistent with typi- 
cal moderate- luminosity AGNs in ULIRGs/SMGs, while 
two of our objects (AzLOCK.Ol and AzLOCK.lO) are 
more consistent with extreme AGN-dominated types. 

3.2.2. Mid-infrared continuum (AGN) properties 

In light of the results in § 3.2.1, we assume from now 
on that the continuum components of the fits are due 
to AGN activity, and we compare the mid-infrared AGN 
properties of our sample with more typical SB-dominated 
SMGs. In general, the continuum portion of our fits from 
§ 2.4 reveals that our SMGs display more strongly rising 
mid-infrared continuum when compared to more typi- 
cal PAH-dominated SMGs (which comprise ~ 80% of 
the SMG population; P08; MD09). In addition, while 
we do not see strong evidence for silicate absorption at 
'^9.7 fim in our IRS spectra (which is often seen in AGN- 
dominated source), it cannot be ruled out since the S/N 
of the continuum between 9-10 ^m is not high enough. 
Notably, all of the best-fitting models to our sample re- 
quire some level of 'additional' continuum component 
compared to typical SMGs. Unfortunately, the data are 
of relatively low S/N and possess insufficient wavelength 
coverage to provide robust measurements of individual 



object spectral slopes. Nevertheless, for completeness, 
we find a range of a = 0-4.6 (uncorrected for extinction, 
and assuming the convention Si, oc with a mean 

and median of 1.6 and 1.8, respectively. In contrast, 
MD09 report a median a of 1.05 for their radio-detected 
sample of (mostly) PAH-dominated SMGs. In addition, 
P08 found that 8/13 of their sample of more typical SB- 
dominated SMGs require a negligible level of continuum 
emission (and more often none) compared to the PAH 
template that was fitted. 

3.2.3. PAH properties 

Here we present the properties of the PAH features 
measured in § 2.4.1 (see Table 4), a key signature of lu- 
minous SF galaxies. We have tried to achieve the best 
PAH line measurements possible for this sample in order 
to compare our results with those of more typical SB- 
dominated samples of SMGs, and to other galaxy pop- 
ulations such as ULIRGs and AGNs, in order to yield 
insight into the physics of these systems. However, many 
of our PAHs are swamped by continuum emission and so 
this makes a fair comparison with other samples difficult, 
given the intrinsic differences in the spectra, the system- 
atic differences in the way lines are measured for the 
different sub-samples (e.g. how extinction is taken into 
account), and the large measurement uncertainties in the 
derived quantities (see e.g. P08 for a discussion). Nev- 
ertheless, qualitatively we find that within the system- 
atic uncertainties our 6.2 fj,m PAH luminosities^ (median 
2.7 X 10^" Lq) are towards the top end of SB-dominated 
SMGs studied by P08 and MD09, who find a spread cov- 
ering about ^ 0.5-3.5X 10^° L0 in their samples. The 
errors on the EWs are currently too uncertain to deter- 
mine if there is a genuine spread in the EW of SB/SMGs 
compared with our sample, although we see the trend 
noted by P08 and MD09, that the SMGs in our sample 
with the largest AGN fractions have the smallest 6.2 and 
7.7 /im EWs. Next, we discuss the similarities and dif- 
ferences in the bulk PAH properties between this sample 
and others by constructing a composite spectrum of our 
AGN-dominated SMGs (see § 3.3). 

3.3. Composite for AGN-dominated SMGs 

The spectral decomposition has revealed that some of 
the SMGs in our sample look more similar to those in 
P08, with low AGN contributions to the mid-infrared, 
while some of our SMGs clearly look more like the rarer 
AGN-dominated Compton-thick objects like CI from 
P08. In order to compare the ensemble of properties of 
AGN-dominated versus SF-dominated SMGs and other 
galaxy populations, we have combined the IRS spectra 
of AGN-dominated (> 50% continuum contribution in 
the mid-infrared) SMGs from all available samples and 
created a higher S/N composite spectrum. Using this 
composite we can search for and characterize faint fea- 
tures. 

The composite is comprised of ten SMGs including 
our five mid-infrared AGN-dominated SMGs, CI (82%) 
and GN04 (61%) from P08, as weU as SMMJ105238.30, 

^ We limit our comparison of PAH line fluxes to the 6.2 /im fea- 
ture since the 7.7 /im feature is blended with other PAH features 
(8.6 /xm, for example) and the 9.7/im extinction trough which fur- 
ther complicates comparisons between different samples. 
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Fig. 7. — Composite IRS spectrum of 10 mid-infrared AGN- 
dominated SMGs (see text). The shaded region shows the Icr 
uncertainty. The red dashed Hne is the best-fitting model (see 
§ 2.4.2), which is composed of a PAH template (green dot-dot-dot- 
dashed line) and a power-law component with extinction applied 
(blue dot-dashed line), revealing a 72% AGN fraction and an un- 
absorbed power-law component slope of ~ 2 (see § 3.3). 
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Fig. 8. — Composite IRS spectrum of 10 mid-infrared AGN- 
dominated SMGs (black histogram, with the shaded region repre- 
senting the 1 cr uncertainty) . We have overlaid the SB-dominated 
SMG composite from P08 (dashed curve), and note that the AGN- 
dominated SMG composite is ~ 5 times more luminous at ^ 7 fim 
than the SB-dominated SMG composite). We have also overplotted 
the T9.7 < 1 weak PAH composite of 17 24 /xm-selected ULIRGs 
(dotted histograms) from Sajina et al. (2007) twice, once renor- 
malised at ~ 7 fim to match the AGN-dominated composite and 
also separately to match the SB-dominated composite, in order 
to highlight the relative similarities and differences between the 3 
composites. This comparison reveals that AGN-dominated SMGs 
appear to have features broadly in between that of SMGs and 
ULIRGs. 

SMMJ123553.26, and SMMJ163650.43 from MD09^ 
Since the individual spectra span a wide dynamic range 
in rest-frame luminosity, we normalized them all to the 
median 7/im luminosity of the sample^. We calculate 
the standard deviation of the sample as a measure of 
the error on the composite spectrum, which dominates 
over the uncertainty from individual flux elements. We 

^ Since MD09 do not explicitly provide individual mid-infrared 
AGN fractions, we have re-derived them here, following the pro- 
cedure in § 2.4.2, and find AGN fractions of 94, 85, and 69%, 
respectively, fulfilling our 'AGN-dominated' criterion. 

''' We normalise the spectra to the median total luminosity be- 
tween rest-frame 6.5-7.5 fim, where there is a lack of PAH features 
or absorption, in order to minimize the potential bias. 



restrict the final SMG composite spectrum wavelength 
coverage to areas where there are > 3 data files, resulting 
in a coverage of '--'4-10 /im in the rest-frame (see Fig. 7). 

We perform a spectral decomposition on the AGN com- 
posite, following the method outlined in § 2.4.2, revealing 
a 72% AGN fraction and an unabsorbed power-law com- 
ponent slope of 2. The amount of extinction required is 
T9.7 ~ 0.5, but is not well constrained (as mentioned pre- 
viously) since the spectrum does not extend much past 
9.7 (but by eye, it seems that the slope is flatter due 
to extinction at longer wavelengths). 

In Fig. 8 we compare our mid-infrared AGN-dominated 
SMG composite to the more typical PAH-dominated 
SMG composite of P08 and to a composite oi z ^ 2 
bright (5*24 > 0.9 mJy) red 24 /im-selected ULIRGs from 
Sajina et al. (2007), which are mainly weak- PAH emit- 
ters and AGN-dominated sources. First of all, Fig. 8 
demonstrates that our AGN-dominated composite is ~ 5 
times brighter at ^ 7 than the SB-dominated SMG 
composite, and we have not renormalizcd them to the 
same luminosity in order to highlight this difference. 
Given that the observed mid-infrared spectra and fluxes 
are so similar between these two samples, about a fac- 
tor of ~ 1.5-2 can be accounted for simply by the ap- 
parently higher average redshift of the AGN-dominated 
combined sample (median z — 2.56) compared with the 
SF-dominated SMG sample (median z — 1.98), and we 
discuss the possibility of selection effects further in § 4. 
Fig. 8 reveals that on average these SMGs have propor- 
tionately stronger/steeper mid-infrared continuum emis- 
sion than more typical SF SMGs, as well as PAH emission 
at 6.2, 7.7, and 8.6 ^m, as well as a hint of silicate ab- 
sorption at ~ 9.7 fim (although the data are very noisy 
in this region of the spectrum). At > 9-10 /xm the av- 
erage AGN-dominated SMG is more similar to the low 
opacity (r9.7 < 1) 24 /im-selected ULIRGs from Sajina 
et al. (2007) than to the SMG composite. In summary, 
in terms of their mid-infrared spectra, AGN-dominated 
SMGs broadly appear to have features somewhere in be- 
tween SMGs (strong SF activity indicated by PAH emis- 
sion) and 24 /im-selected ULIRGs (ongoing SF with a 
relatively strong AGN component). 

3.4. Full SED Fits 

We have established mid-infrared AGN-dominance for 
the majority of our SMG sample. However, it is not 
clear if the AGN activity will continue to produce a 
dominant fraction of the total infrared luminosity, since 
it is expected that the AGN will become less impor- 
tant bolometrically than SBs when extrapolating from 
the mid-infrared to longer wavelengths (Tran et al. 2001; 
P08). Here we investigate the effect of the AGN emission 
to the total infrared luminosity when longer wavelength 
(submm) emission near the peak of the bolometric energy 
output of the galaxy is taken into account. 

To begin with, we fit the mid-infrared spectrum and 
submm/mm data (but not the radio data; see § 3.4.1) 
for each source with the full SF-SMG composite SED 
from P08, which describes typical SMGs well (and in- 
cludes a ~ 30% AGN contribution). Visually, all of the 
fits look poor and unsurprisingly have values of x^cd ^ 1' 
meaning that the assumed model does not provide a very 
good description of the data. Since we already have evi- 
dence from the mid-infrared fitting in § 3.2.1 supporting 
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Fig. 9. — Mid-infrared-to-submm SEDs of our sample of SMGs. The dark solid curve shows the best fitting linear combination of an 
SB component (the SMG starburst composite SED from P08 comprised of their SMG SB composite spectrum spliced with the best-fitting 
modified CEOl template, dotted curve) and an AGN component (a Mrk231 template, dashed curve) to the IRS spectrum (histogram) 
and 850 or 1100 /.im photometry. Source names are given in each panel (see Table 1), as well as the best-fit values. We also include 
350 /.tm photometry (Coppin et al. 2008a) in the fitting for LOGK850.01. We note a significant degradation in the quality of the fits in the 
mid-infrared if the additional AGN component is not included, indicating that this subset of SMGs is generally well described by the SEDs 
of typical SMGs, but with an additional AGN component required. For completeness, we overlay 3.6— 24 /^m and 610 MHz and 1.4 GHz 
photometry (see Table 2) as well as 70 and 450 /.im 3cr upper limits from Hainline et al. (2009) and Goppin et al. (2006), respectively, to 
demonstrate that these data are roughly consistent with the template fitting, even though they are not used explicitly. 



a significant continuum/ AGN contribution over typical 
SMGs in many cases, we repeat the fitting with a linear 
combination of the SF-SMG composite and an additional 
AGN component represented by Mrk 231, as described in 
§ 2.4.3. Visually the fits improve drastically, and we find 
an improvement in Xrod values in every case, typically of 
Axred -3 (though as large as ~ 10 for AzLOCK.lO and 
^ 60 for AzLOCK.Ol), supporting the inclusion of the 
additional AGN template into the total SED model (see 
Fig. 9). 

A wider range of templates and/or models could be 
explored in principle, however many complications could 
arise and make the results difficult to interpret. Since we 
wish to compare our decomposition directly with that of 
more typical SB-dominated SMGs (P08), the most natu- 
ral choice for the AGN component is the SED of Mrk 231 
since this is also what P08 used. Mrk 231 is a well-studied 
local ULIRG with strong evidence for an emerging AGN. 
For example, Genzcl et al. (1998) conclude from high- 
excitation line measurements that Mrk 231 is likely to be 



powered mainly by an AGN, and more recently, Armus 
et al. (2007) find a mid-infrared AGN fraction of - 100%. 
Although Mrk 231 appears to be AGN-dominated in sev- 
eral wavebands it could also harbor some amount of 
dusty SF which would bias our bolometric AGN fractions 
upward. In order to verify this possibility we repeat the 
SED fitting with the nuclear spectrum of the well-known 
local AGN NGC 1068 (R. Chary priv. comm.; Genzel et 
al. 1998; Lutz et al. 2000) in place of Mrk 231 and, al- 
though the relative shapes of their dust spectra indicate 
that NGC 1068 is cooler than Mrk 231, the fitting yields 
the same AGN fractions to within a few percent. We 
also explore SF / AGN combinations using M82 as the SB 
component (instead of the P08 composite) and find that 
while it reproduces the mid-infrared PAH features well, 
M82 catastrophically fails to fit the submm/mm pho- 
tometry, as expected since M82 has a much warmer dust 
temperature than SMGs. A more extensive comparison 
using a suite of many-parameter models (e.g. Siebenmor- 
gen & Kriigel 2007) is beyond the scope of our simple ap- 
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preach given the sparsity of the data, ahhough could be 
attempted in the near future as the data quality and SED 
coverage could be improved with upcoming Herschel snd 
JWST surveys, for example. 

Using these fits, we have calculated values of total Lib, 
as well as Lsb and -Lagn for the SF and AGN compo- 
nents separately and give these in Table 4 (sec § 2.4.3). 
The fractional contribution of AGN to the total Lir in 
our sample ranges from 17-63%, with evidence for only 
two SMGs being bolometrically dominated by an AGN 
(> 50% maximal AGN contribution - AzLOCK.Ol and 
AzLOCK.lO which also have the strongest AGN com- 
ponents in the mid-infrared). Combining our five SMGs 
with a > 50% AGN-contribution with CI and GN04 
from P08, the average AGN contribution to the bolo- 
metric luminosity, Lir , is ~ 40%, compared with 15% for 
typical star-formation dominated SMGs from P08. We 
thus find a median Lsb for our sample of 6.5 x 10^^ Lq, 
which translates into a median SFR for the sample of 
~ 1000M©yr~^, similar to typical bright SB-dominated 
SMGs with 5 x IO^^l^ (P08). Our sample thus has a 
similar Lm to typical SB-dominated SMGs (only 1.5 
times more compared to P08), and thus it appears that 
our SMGs are from a similar SB luminosity class with 
similar SFRs, but that our sample has proportionately 
higher AGN content in the mid-infrared. We discuss the 
implications of these results and how they might fit into 
an evolutionary context with other SMGs in § 4. 

We have only fit the SEDs to the mid-infrared spectra 
and the submm data, but note that extrapolation of the 
best-fitting SEDs does generally fit the observed IRAC 
and 610 MHz photometry of our SMGs well, with some 
amount of scatter at 1.4 GHz (see Fig. 9). Radio fre- 
quencies of 1.4 GHz arc sensitive to synchrotron radiation 
from relativistic electrons from supernovae (and hence 
recent star-formation; Condon et al. 1992), although an 
AGN can also produce comparable synchrotron emission 
by jets and radio lobes, leading to ambiguities in the 
cause of the mechanism at /iJy levels in the absence of 
other information. Since we know that these SMGs are 
harboring AGN, it is quite possible (and likely) that the 
cause for the disagreement between the SED fits and the 
1.4 GHz data is due to an AGN which we now investigate 
by examining their measured 1.4GHz-610MHz spectral 
indices for clues. 

3.4.1. Radio Spectral Indices 

Ibar et al. (2010) calculate 1.4GHz-610MHz radio 
spectral indices based on radio maps from Ibar et al. 
(2009) using a flxed-to-beamsize source extraction, and 
they find a mean 1.4GHz-610MHz radio spectral index 
of ali% = -0.75±0.06 (5*^ oc i^") for the bulk of SMGs in 
the LH, consistent with optically thin synchrotron emis- 
sion. The spectral indices for those of our SMGs which 
lie in the LH have a range of —1.32 < a < —0.05 (see 
Table 2). Five of these LH SMGs have steep radio spec- 
tra, al'iQ <, —0.84, while one has a self-absorbed flat 
spectrum, all of which arc clearly deviant from typical 
SMGs. This difference seen in this subset of SMGs is 
suggestive that the radio emission has a diflFerent origin 
in these systems. Since we know from the mid-infrared 
spectral signatures that these SMGs likely harbor AGN, 
it seems reasonable to assume that the AGN is affecting 
the radio emission, although other environmental causes 



could be playing a role (Ibar et al. 2010). In order to 
quantify this we calculate the logarithmic ratio of the 
rest-frame infrared and radio luminosities for SMGs us- 
ing the parameter from Kovacs et al. (2006), and find 
a range of values from 2.0-2.5, with a median value 
of 2.21 ± 0.15 compared to 2.14 ± 0.07 for radio-detected 
bright SMGs. Although there is no evidence for an off- 
set in <7l compared to Kovacs et al. (2006), the spectral 
indices suggest that the radio emission could be contam- 
inated by radio-quiet AGN emission, providing some jus- 
tification for omitting these data from the fitting. 

4. DISCUSSION 

These new IRS spectra of AGN-dominated SMGs have 
helped to confirm the predictive power of the S24/SS.0- 
5'8.o/'S'4.5 diagram. Our sample was selected to lie off 
the expected SB sequence in the diagnostic color-color 
plot (Fig. 1), and we have confirmed that the contin- 
uum emission dominates over PAH emission in the mid- 
infrared for 5/8 of our sources, which we interpret as due 
to a powerful AGN. The majority of our SMGs possess 
stronger/steeper mid-infrared continuum emission than 
more typical SB-dominated SMGs, indicating that there 
is a continuous trend in AGN content in SMGs - from 
low-luminosity AGN to extremely luminous AGN. As 
with all color-selection methods, the IRAC color selec- 
tion is merely a guide to singling out the most powerful 
AGN, and some level of contamination from SB systems 
is to be expected. In particular, SF systems at 2; < 1 
could have 8^/84,5 > 2, since the 4.5 and 8^m chan- 
nels are not yet climbing up the restframe 1.6 /xm stellar 
bump. For example, the most discrepant data point on 
the AGN track is an SMC at z = 0.689 from Hainline 
et al. (2009) which is clearly SB-dominated in the IRS 
spectrum (MD09). The color cut is thus appropriate for 
separating mid-infrared SB- and AGN-dominated SMGs 
from z ~ 1-4. However, since most SMGs lie within 
1 < z < 4 (Chapman et al. 2005), negligible contamina- 
tion is expected on either side of the color cut boundary. 
In light of the IRS results of MD09, who observed a large 
fraction of the submm-bright radio-selected sample of 
Chapman et al. (2005) which have optical spectroscopic 
redshifts, we advocate adjusting the selection slightly 
to Ss/S4,5 > 1.65 to account for the AGN-dominated 
SMGs seen in the MD09 sample. Even then, we still see 
some small amount of scatter to either side of this line 
(see Fig. 1). Combining all the IRS-observed blank-field 
SMC samples, we find that ~ 15% of SMGs are mid- 
infrared AGN-dominated systems. Our full SED fitting 
reveals that 2/8 (25%) of SMGs in our sample with AGN- 
dominant mid-infrared emission are likely bolometrically 
dominated by an AGN (in the far-infrared). This im- 
plies that overall ~ 5% of blank-field SMGs are likely to 
have a dominant contribution to their Lir due to AGN 
activity. 

Now that we have confirmed that the majority of our 

sample is comprised of star-forming SMGs with an ex- 
cess of AGN emission compared to typical SMGs, we 
can interpret this result within the framework of the 
proposed evolutionary sequence which links our AGN- 
dominated sample of SMGs to typical SMGs. However, 
before placing these results in any wider context, we must 
first take stock of any potential selection effects and how 
they could bias our interpretation. 
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4.1. Selection Effects 

The eight targets in our color-selected sample have 
850 /im and 24/im flux densities of 4-13 mJy^ and 180- 
1600 n3y, with median values of 9 mJy and 310 /iJy, re- 
spectively (see Table 2). Note that our targets have 
a similar 24/im flux density distribution as the whole 
SHADES SMG sample (i.e. our color selection has not 
biased the sample to preferentially brighter or fainter 
24 /xm SMG counterparts). Thus, our IRAC color- 
selected subset of SMGs studied here should be repre- 
sentative of ~ 15% of the SMG population discovered 
in existing blank-field 850 /itm imaging. The primary 
submm/mm-selection uniformly identifies high-redshift 
galaxies with Lsb > 10^^ Lq and corresponding SFRs of 
~ lOOOMeyr^i. The AGN-dominated SMG sample we 
have focussed on here is similar to typical SMGs in terms 
of Lsb , which allows for a fair and direct comparison be- 
tween the two sub-samples of SMGs. However, there is 
some evidence that our sample lies at a slightly higher 
average redshift than typical SMGs, which merits some 
discussion. 

The redshift distribution of our color-selected sample 
of eight objects ranges from 2.5-3.4, with a median of 2.8. 
This is noticeably shifted to higher redshifts compared to 
the P08 and MD09 samples, which have a median red- 
shift of 2.0, more consistent with the known redshift dis- 
tribution for radio-identified SMGs with Ss5o > 5mJy 
{z ^ 2.2; Chapman et al. 2005). The bias in our red- 
shift distribution to earlier epochs is likely a submm 
selection effect, which tends to pick out warmer galax- 
ies on average at higher redshift (e.g. Pope et al. 2006; 
Coppin et al. 2008a; Dye et al. 2009). The subsequent 
IRAC color-selection imposed could be eflfectively sin- 
gling out those higher-redshift SMGs which possess warm 
dusty tori heated by a proportionately larger AGN. This 
trend was also noted by MD09 and Watabe et al. (2009), 
who found that 'typical' SMGs show an increase in mid- 
infrared continuum strength (and less prominent PAHs) 
with increasing redshift. One way to verify this trend is 
to measure the flux near the far-infrared dust peak of the 
SED in order to constrain the dust temperatures in our 
Ss/Si.s > 2 sample and compare them to the dust tem- 
peratures of more typical SMGs with Ss/S4.5 < 2, over 
a wide redshift range. Facilities such as Herschel and 
SCUBA-2 will be able to help address this issue for this 
sample, although some progress could be made now using 
BLAST which has already measured the rest-frame IR 
peaks of hundreds of SMGs in other fields where Spitzer 
IRAC data exist (Devlin et al. 2009; Dye et al. 2009). 
However, it is unlikely that the average redshift differ- 
ence (equivalent to ~ 1 Gyr in cosmic time) will have a 
large effect on the main conclusions, since even SMGs at 
z ~ 4-5 seem to show a mix of starburst and AGN signa- 
tures/SB properties similar to SMGs at 2; ~ 2 (e.g. Cop- 
pin et al. 2009). 

4.2. Where do our AGN-dominated SMGs fit into the 
evolutionary sequence? 

* Note that wc have converted the AzTEC 1.1mm fluxes to 
850 ^m fluxes by multiplying by a factor of 1.8, according to the 
average 5850/51100 color ratio observed for SMGs (e.g. Chapin et 
al. 2009). 



For the purposes of this discussion, we adopt the pro- 
posed evolutionary sequence first proposed by Sanders et 
al. (1988) and test the scenario where the SMG popula- 
tion evolves through a primarily starburst phase, onto a 
subsequent submm-detected AGN-phase, into a submm- 
undetected QSO and finally a passive elliptical. Spitzer- 
IRS data has tentatively identified some local ULIRGs 
being in key stages in this lifecycle (Farrah et al. 2009b) . 
Now we ask - are the observed properties of typical SMGs 
and the AGN-dominated SMG sample studied here con- 
sistent with this picture in the high-redshift universe, in 
terms of their space densities and other observed proper- 
ties? 

Ss5o > 4.5 mJy SMGs have surface densities of ~ 

600 dcg~^ (Coppin et al. 2006) or space densities of 
~ lO"-"^ Mpc"'^ for z = 1 3 (Chapman et al. 2005; Wall, 
Pope & Scott 2008). The combination of this IRS work 
and other samples of SMGs has uncovered a mid-infrared 
presence of appreciable heating of a dust torus by an 
AGN for approximately 15% of the SMG population at 
this epoch. Assuming that all bright SMGs go through a 
subsequent AGN phase then the observed IRS-confirmed 
AGN-detection fraction merely represents a 'duty cy- 
cle'. The current total lifetime estimates from modelling 
SMGs are ~ 100-300 Myr (Swinbank et al. 2006), and 
thus the inferred submm-bright AGN-dominated phase 
caught by the mid-infrared would be occuring during 
the last ~ 15-45 Myr of an SMG hfetime. These rel- 
ative lifetimes and duty cycle are consistent with the 
merger simulations of Springel, Di Matteo & Hernquist 
(2005), where the tidal forces from galaxy mergers trig- 
ger a nuclear SB, fuelling a rapid growth of BHs before 
the gas has been consumed. The BH accretion rate grows 
throughout the peak phase of SF in the merger system 
(with the peak phase lasting « 100 Myr in these simula- 
tions) , peaking « 50 Myr later than the SFR peak due to 
the delayed action of AGN feedback on the gas and man- 
ifesting itself as a luminous QSO for w 50 Myr total. In 
the Springel, Di Matteo & Hernquist (2005) simulations, 
the SMG phase is in principle both an SB and an AGN, 
but it is only during the final evolutionary stage that 
the remnant would be visible as an AGN, with outflows 
removing the surrounding gas and dust which heavily ob- 
scures the AGN at the beginning of the burst. Tentative 
evidence for this has been found in typical SMGs since 
their estimated MBH/-M^gai ratio is lower than found in 
z^2 QSOs (Alexander et al. 2008a; Coppin et al. 2008b; 
Peng et al. 2006) 

Recently, Dey et al. (2008) proposed that Dust Ob- 
scured Galaxies (DOGs; S^{24iim)/S^{R) > 1000 and 
^24 > 300 mJy, with surface densities similar to > 6 mJy 
SMGs) could represent a brief evohitionary phase be- 
tween SMGs and less obscured QSOs or galaxies, which 
is supported by the finding that ^ 30% of SMGs (includ- 
ing the AGN-dominated SMGs) meet the DOG criteria 
(Pope et al. 2008b). Interestingly, we note that virtually 
all (5/7 with i?-band coverage) of our targets here fulfill 
the DOG criteria, confirming that the majority of AGN- 
dominated SMGs also overlap with the DOG population. 

While these arguments are mainly circumstantial, we 
have shown that the AGN-dominated SMGs at least fit 
into the evolutionary picture of Sanders et al. (1988) 
in terms of the inferred duty cycles. While the sta- 
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tistical approach to duty cycles and relative lifetimes 
can be useful, in-depth studies of supposedly individual 
SMGs often reveal a complicated picture involving mul- 
tiple components (e.g. Ivison et al. 2008). Measuring gas 
masses with the IRAM Plateau de Bure Interferometer 
or the Atacama Large Millimeter Array, inferring gas de- 
pletion timescales, and estimating SMBH masses using 
near-infrared spectroscopy for the AGN-dominated sam- 
ple studied here (see e.g. Alexander et al. 2008a; Coppin 
et al. 2008b) would allow us to put further constraints 
on the SMG evolutionary sequence we are exploring. 

5. CONCLUSIONS 

Spitzer IRS spectroscopy has been obtained for a sam- 
ple of eight SMGs brighter than 200 /iJy at 24 /im and 
selected to possess observed IRAC colors of 5'8/S'4.5 > 2, 
in order to investigate the level of AGN dominance in 
the SMG population. Our conclusions are as follows: 

1) Taking into account sources in the literature, we 
find a that a color-selection of Ss/S4,5 > 1.65 is a bet- 
ter description overall for defining the boundary between 
SB and AGN-dominated SMGs, with a small amount of 
scatter across this division. 

2) We see signs of both SF and AGN activity in our 
sample of SMGs, with a continuous distribution of AGN 
fractions in the mid-infrared. There are signs of SF in the 
form of PAH features in all of the SMGs in oiu- sample, 
from which we derive redshifts between 2.5-3.4, demon- 
strating the power of the mid-infrared to determine red- 
shifts when the optical counterparts are too faint to study 
with current facilities. 

3) We find a median AGN contribution to the mid- 
infrared of 56% in OUT sample of SMGs, indicating that 
overall, SMGs with > 0.2 mJy and 58/5*4.5 > 2 tend 
to have more dominant AGN-components in the mid- 
infrared than typical SMGs (< 30% mid-infrared AGN 
contribution). Extrapolation to the far-infrared reveals 
that the AGN is bolometrically important only in two of 
our SMGs (each with a > 75% mid-infrared AGN con- 
tribution). This result suggests that significant contami- 
nation from AGN to the far-infrared luminosities affects 
g 5% of the SMG population overall. 

4) To first order, these AGN range from being fairly 
low-luminosity types found in typical SMGs to more 



extreme cases that likely overlap with AGN-dominated 
24 /im-selected ULIRGs in Sajina et al. (2007). In any 
case, all our AGN appear to be obscured to some degree 
(at least five have low ix/ie/im and six are optically 
faint with R > 23) and a small subset of our sample may 
be obscured by Compton-thick material. 

5) Our results are thus consistent with the Sanders et 
al. (1988) evolutionary scenario, with all SMGs under- 
going a 'transitional' AGN-dominated phase with a duty 
cycle of ~ 15%. Our sample of AGN-dominated SMGs 
could be at a slightly later stage of evolution than SF- 
dominated SMG systems, with the SF still occurring but 
where the AGN has now begun to heat the dust appre- 
ciably in the SMG as the BH undergoes a period of rapid 
growth. 
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TABLE 1 

Spitzer-IRS OBSERVATIONS OF SHADES SMGs in the Lockman Hole ('LOCK') and the Subaru- MM Deep 

Field f'SXDF'). 
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Nickname 
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Integration time 








R.A. 


Dec. 


LL2(xl20s) 


LLl (xl20s) 


Total (hr) 
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2x6 
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10x6 


4.8 
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20x6 


16x6 
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TABLE 2 

Optical-to-radio photometry of the SMG 24 ^6M counterparts 



Nickname 
AzLOCK.Ol 


B 


R 




z 


22. 


J 

.84 ± 0.10 


21. 


K 

12 ± 0.04 


3.6 i^m 
42.1 ± 4.5 


4.5 nm 
62.1 ± 5.8 


2362.2 ± 79.9 


8 fiui 

(f Jy) 

194.8 ± 14. 


3 


16 i^ia 

(MJy) 

697.5 ± 51.5 


24 (im 
(AiJy) 
1587.3 ± 43.6 


850 or 1100 Mm 

(mJy) 
eoli.o 


, 1.4GHz 

(cJy) 

258 ± 11 


„350 
"1.4 

0.69 ± 0.03 


^1.4 
"610 

-1.01 ± 0.07 


SXDP850.02 


26.12 


25.33 


25.26 


24.96 




24.40 




22.54 


17.0 ± 0.8 


22.5 ± 0.8 


< 48 


60.9 ± 4. 


5 


< 171.0 


313.0 ± 47.0 


10 2+1-^ 
^"■''-1.6 


66 ± 11 


0.91 ± 0.06 




AzLOCK.lO 


23.04 ± 0.01 


22.65 ± 0.01 


22.40 ± 0.01 


22.23 ± 0.01 


20 


.87 ± 0.06 


20. 


16 ± 0.02 


39.7 ± 4.4 


50.1 ± 5.2 


71.0 ± 14.0 


116.6 ± 11. 


1 


269.1 ± 46.1 


668.4 ± 29.5 


4 1 + 0.9 


77 ± 9 


0.83 ± 0.06 


-1.32 ± 0.20 


LOCK850.15 


26.54 ± 0.12 


25.02 ± 0.05 


25.24 ± 0.09 


24.50 ±0.11 




> 22.33 


22. 


25 ± 0.12 


9.7 ± 2.2 


15.4 ± 2.9 


< 18 


36.7 ± 6. 


5 


105.4 ± 46.0 


353.0 ± 20.0 


^^■^-5^0 


105 ± 6 


0.88 ± 0.09 


-1.10 ± 0.10 


SXDF850.37 


27.82 ± 0.26 


26.87 ± 0.16 


26.96 ± 0.20 


26.49 ± 0.32 


25 


20 ± 0.79 


23 


25 ± 0.21 


14.8 ± 0.6 


25.2 ± 0.9 


< 48 


75.0 ± 5. 


2 


138.0 ± 53.8 


183.0 ± 47.0 


-«i^2:?. 


41 ± 9 


0.85 ± 0.15 




AzLOCK.05 


25.94 ± 0.07 


25.46 ± 0.08 


25.85 ± 0.16 


25.18 ± 0.21 




> 22.33 




> 22.93 


8.9 ±2.1 


12.5 ± 2.6 


< 18 


29.7 ± 5. 


9 


< 121.8 


235.0 ± 24.1 


-«ti:°9 


138 ± 9 


0.75 ± 0.05 


-0.84 ± 0.11 


LOCK850.01 


> 27.17 


25.83 ± 0.11 


25.35 ± 0.10 


24.89 ± 0.16 




> 22.33 




> 22.93 


5.1 ± 1.6 


8.4 ± 2.2 


12.2 ± 2.2 


17.3 ± 4. 


7 


< 125.7 


217.0 ± 16.0 




110 ± 6 


0.79 ± 0.03 


-0.05 ± 0.18 



Optical— to— radio mult i wavelength magnitudes or flux densities of our sample of SMG 24jum counterpart targets, listed in order of redshift (see Table 3). Non-detections are indicated by the 
corresponding 6(7 point source sensitivity limit of the appropriate data set. Columns 2—5 list the Subaru BRi' z optical photometry (Dye et al. 2008; Furusawa et al. 2008; Mortier et al. in preparation), 
while columns 6—7 provide the Ji^-band photomcrtry from the DR3 UKIDSS release (Lawrence et al. 2007 and Warren et al. in preparation; see § 2.1). The optical/near-IR photometry is in AB 
magnitudes and has been measured in 2" apertures (with 5 cr upper limits given when the extracted photometry is fainter than these limits). The 3.6-8 /^m photometry for the LH SHADES, LH 
AzTEC, and SXDF SHADES sources are from Dye et al. (2008) and Hainline et al. (2009), Egami et al. (priv. comm.), and SWIRE Londsdale et al. (2004) and SpUDS (Dunlop et al. priv comm.), 
respectively, and arc listed in columns 8—11. New 16 /xm photometry data (this work) are listed in column 12. The dcboosted submm and mm photometry from Coppin et al. (2006) and Austermann 
et al. (2010), respectively, are listed in column 14. Columns 13 and 15 give the 24 fim and radio flux densities from Ivison et al. (2007), Ibar et al. (2009, 2010), and Arumugam et al. (in preparation). 
The submm-to-radio spectral indices (350 GHz-to-1.4 GHz) are listed in column 16, with Ss50 for the 1100 /J.m AzTEC sources calculated assuming an 5'85o/5'iioo ratio of 1.8 (see Chapin ct al. 2009). 
The 610 MHz-to-1.4 GHz spectral indices from Ibar et al. (2010) are given in column 17 and are available for the LH sources only. '— ' indicates no coverage at the source position in this band. 
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TABLE 3 

Mid-infrared and X-ray derived properties. 





^IRS Continuum 
[%] 


[10*6 erg s-i] 


[10 


(0 5-2 kpV~l 

ergcm"^ s~^] 


J-^ \^ J-'-' JVC V 1 

[10*3ergs-l] 


AzLOCK.62 


2.48 ± 0.03 


32 


0.416 




< 5.23 


< 3.55 


AzLOCK.Ol 


2.50 ± 0.02 


71 


6.76 (7.21) 




< 9.07 


< 6.31 


SXDF850.02 


2.55 ±0.06 


10 


0.0374 




< 10.95 


< 7.81 


AzLOCK.lO 


2.56 ± 0.03 


77 


2.96 




4.40 


3.19 


LOCK850.15 


2.76 ±0.01 


62 


1.94 




< 6.08 


< 5.16 


SXDF850.37 


2.78 ± 0.04 


51 


0.641 




< 15.14 


< 13.0 


AzLOCK.05 


2.82 ± 0.02 


26 


0.378 (1.84) 




< 6.94 


< 6.15 


LOCK850.01 


3.38 ± 0.02 


56 


1.29 




< 4.78 


< 5.94 



In column 2 we list the spectroscopic redshifts derived from the fitting procedure described in § 2.4.1 
for the PAH features in the mid-infrared spectra. In column 3 we give the AGN contribution to the mid- 
infrared, expressed as a percentage of the total mid-infrared luminosity covered by the IRS spectra, as 
derived from the continuum power-law fits from § 2.4.2. Column 4 gives the resl-framc 6 /im luminosities 
based on the best-fitting power-law component, assuming that the AGN contributes to all of the emission 
at 6 /J.m (with the extinction-corrected values given in parentheses where appropriate). Column 5 gives 
the observed-frame X-ray fluxes of the detection and 3 a upper limits for the non-detections, and column 
6 gives the corresponding rest-frame X-ray luminosities, calculated using the X-ray flux and the IRS- 
derived redshifts; V = 1.4 is used to make small K-corrections (see § 2.1). 
^The counterpart is l" away from the 24jum position. 



TABLE 4 

PAH Luminosities, Equivalent Widths and Infrared Luminosities 



SMO ID 


PAH 


huiiinosity (lO" 






PAll (Xiiiirakait width 


" (inn) 












().2 //111 


7.7 //111 


8.() //111 




().2 //111 


7.7 //111 




8.( 


) //111 


SI! . 




l/jtal 


AzLOCK.62 


19.4 ± 1.1 


82.7 ±2.3 


16.8 ±4.6 


0.62 


±2.30 


3.43 ±11.44 


0.77 


± 


2.48 


3.2 


1.8 


5.0 


AzLOCK.Ol 


76.5 ± 24.5 


387.1 ± 11.1 


64.4 ±5.1 


0.26 


±0.22 


1.80 ±1.68 


0.34 


± 


0.33 


7.9 


13.6 


21.5 


SXDF850.02 


15.8 ± 9.4 


161.9 ±8.0 


22.9 ±11.5 


0.62 


±5.03 


6.61 ± 37.52 


0.96 


± 


4.50 


8.3 


1.6 


10.0 


AzLOCK.lO 


30.1 ± 14.9 


168.4 ± 28.7 


26.4 ± 11.2 


0.23 


±0.14 


1.59 ±0.84 


0.27 


± 


0.14 


5.5 


5.9 


11.3 


LOCK850.15 


26.6 ± 6.5 


125.2 ±4.2 


36.3 ±4.5 


0.26 


±0.17 


1.23 ±0.52 


0.36 


± 


0.13 


10.7 


4.2 


14.9 


SXDF850.37 


12.6 ±8.1 


86.8 ±4.3 


8.9 ±8.5 


0.42 


±0.88 


4.63 ± 12.19 


0.54 


± 


1.52 


3.6 


1.6 


5.2 


AzLOCK.05 


36.5 ±8.7 


193.0 ± 8.0 


2.8 ±4.5 


1.21 


±8.38 


6.58 ±31.27 


0.10 


± 


0.37 


7.8 


2.3 


10.1 


LOCK850.01 


25.7 ±5.0 


96.9 ±31.9 




0.27 


±0.50 


0.59 ± 0.44 








6.5 


3.3 


9.8 



^Thc equivalent widths are all in the rest-frame. 

^Obtained from the best fitting linear combination of an SB component (the SMG starburst composite SED from P08 comprised 
of their SMG SB composite spectrum spliced with the best-fitting modified CEOl template) and an AGN component (a Mrk231 
template) to the IRS spectrum and 850 or 1100 jLtm photometry (see text). The errors are typically ^ 10% on each component. 



APPENDIX 

NOTES ON INDIVIDUAL SOURCES 

LOCK 8 5 0.1 

LOCK850.01 is the brightest SCUBA source detected in the LH region of the 8mJy Survey (Scott et ah 2002; 
LE850.01) and SHADES (Coppin et al. 2006) and is also one of the brightest 1200 /xm sources detected in the Greve 
et al. (2004) MAMBO survey (LH1200.005). Lutz et al. (2001) identify extended K-hand emission with the PdBI 
1.2mm-confirmed position, and using the combination of redshift estimates based on the photometric properties of 
the Extremely Red Object (ERO) counterpart they provide a 'best guess' of z = 3. Further spectroscopic follow-up 
of LOCK850.01 was obtained through optical and near-infrared spectroscopy (Simpson et al. 2004; Chapman et al. 
2003; Chapman et al. 2005; Blain et al. 2004; Ivison et al. 2005), although no secure redshift was forthcoming from 
the optically faint ERO. A galaxy ^ 3" away from the radio counterpart was also targeted and estimated to lie at 
z = 2.148 based on absorption lines in the LRIS spectrum. This redshift was tentatively assigned to LOCK850.01, 
since the likelihood of finding a z ^ 2 galaxy so close to the SMG ccntroid by chance is slim and so was deemed likely 
to be associated with the optically- faint ERO (sec Ivison et al. 2005 for a complete discussion). Based on different 
photometric redshift techniques Aretxaga et al. (2007) and Dye et al. (2008) find Zphot = 2.4^ 2 (with a 90 per cent 
confidence range of 2.2-3.8) and Zphot = 4.21 (with a 90 per cent confidence range of 2.41-4.45). respectively. Based 
on the presence of two PAHs in Fig. A, z — 3.38 ± 0.02 is the best-fit solution and we adopt this redshift for 
LOCK850.01. 
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LOCK850.41-2 

LOCK850.41 was first discovered in the SmJy Survey (Scott et al. 2002; LE850.08) and by the Greve et al. (2004) 
1200 /xm MAMBO survey (LH1200.014), and was reconfirmed by SHADES with an 850 /xm flux of 3.9t'j':^mJy (de- 
boosted) and a submm position of RA=10:51:59.86, Dec=57:24:23.6 (Coppin et al. 2006). Ivison et al. (2005, 2007) 
and Chapman et al. (2005) located a statistically robust double counterpart in the radio for this system, as well as two 
coincidental 24 /im counterparts, and obtained a convincing redshift of z = 0.689 for the brightest, but furthest, of the 
counterparts. A spectrum for the other radio counterpart was also obtained, containing broadly similar spectral fea- 
tures but at much lower S/N and so a redshift could not be unambiguously determined. The redshift of the z — 0.689 
counterpart (LOCK850.41-1) was confirmed by MD09 using the Spitzer-IRS. To complete the picture we obtained 
an IRS spectrum of LOCK850.41-2. This source originally followed our IRAC color selection criteria, although with 
a proper deblending of the IRAC photometry, it was found not to make our color selection cut after all, and so we 
have excluded this source from the main results and discussion of the paper. In Fig. Al, the source appears to be 
completely continuum-dominated, with no strong discernable PAH features from which to determine a redshift. We 
find a formal 'minimum' best-fit at z — 3.1 (see § 2.4.1) if the weak feature seen at ~ ll/j,m is assumed to be 
the 6.2 /xm PAH feature, although we note that this solution is not unique as all the reduced values are < 1 for 
the 0.5 < z < 4.5 range that was considered. This counterpart possesses a hard X-ray identification ^ 2.9" away 
('S'x (0.5-2 koV) — 17.5 X 10^^^ ergcm"^ s~^; ix(2-iokcV) — 18.3 x lO^'^ergs"^; see Fig. 2) which is spectroscopically 
classified as a Type 1 AGN in the optical, which could explain the lack of visible strong PAH features in the IRS 
spectrum at z = 0.974 (Lehmann et al. 2001). Mainieri et al. (2002) find no evidence for X-ray absorption in this 
source, but do find a surprisingly flat X-ray spectral slope of F ~ 0.7, compared with a typical F ~ 1.9 for Type 1 AGN. 
Despite the potential redshift ambiguity for FOCK850.41-2, the source is continuum-dominated in the mid-infrared, 
regardless of the redshift. 




Rest wavelength (|im) 

Fig. Al. — The IRS spectrum of LOCK850.41-2: a continuum-dominated source in the mid-infrared (< 90% continuum/AGN fraction) 
with no visible PAH features from which a redshift could be obtained. Instead we assume a source redshift of 2: = 0.974, as indicated by 
the optical spectrum of the X-ray counterpart of the SMG from Lehmann et al. (2001). The solid black histogram is the raw (unsmoothed) 
IRS spectrum, while the shaded region represents the associated la noise from the sky background (see § 2.3). The red dashed line is the 
best-fitting model composed of an M82 PAH template (green dot-dot-dot-dashed line) and a power-law component with extinction applied 
(blue dot-dashed Une). See § 2.4.2. 



